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Abstract The past decade has presented a revolution in the field of observational 
high energy y-ray astrophysics with the advent of a new generation in ground-based 
TeV telescopes and subsequent GeV space telescopes. The Fermi Large Area Tele- 
scope (LAT) was launched in August 2008 and has offered unprecedented sensitivity 
and survey capabilities in the 30 MeV - 300 GeV energy range. 

Presented here are the results from the first two years of LAT observations of 
galactic binary systems including the definitive detections of LS I +61°303, LS 
5039 and Cyg X-3. These sources and others are discussed in context with their 
known TeV and X-ray properties. The LAT data provides new understandings and 
pose new questions about the nature of these objects. The identification of an expo- 
nential cutoff in the spectra of both LS I +61°303 and LS 5039 was unexpected and 
poses challenges for explaining the emission mechanisms and processes which are 
in operation within these systems. 
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1 Introduction 

To date there are only five X-ray binaries that have reported detections of high en- 
ergy (HE; 0.1-100 GeV) or very high energy (VHE; >100 GeV) y-ray emission: 
LS I +61°303 [13, 1]; LS 5039 [8, 2]; PSR B1259-63 [9]; Cyg X-3 [15, 3]; Cyg 
X-l [12, 34]. Whilst there were theoretical predictions that these systems may be 
capable of emitting such high energy radiation and some reports of high energy 
sources that may be associated with some of the above systems the first definitive 
detections and identifications have only occurred within the past five years due to 
the current generation of sensitive HE and VHE telescope facilities. 

PSR B 1259—63 was long known to be the only millisecond pulsar in orbit around 
a main-sequence star. It was a known hard X-ray source and was predicted to pro- 
duce y-rays through the interaction of the pulsar wind with the disk of the Be com- 
panion; VHE radiation was discovered by the H.E.S.S. experiment [9] confirming 
this theoretical prediction. LS I +61° 303 and LS 5039 are both persistent, well 
known X-ray binaries which were spatially consistent with HE sources detected 
by the Energetic Gamma-Ray Experiment (EGRET) in the 1990s [25], however the 
lack of detection of periodic flux modulation and the size of the error box meant it 
was impossible to definitively claim a detection of y-rays from these sources. This 
all changed with first the detection of these sources at VHE by the MAGIC, VER- 
ITAS and H.E.S.S. Cherenkov radiation telescopes and subsequently the detections 
at HE by the Fermi Gamma-Ray Space Telescope. Cyg X-l and Cyg X-3 are ex- 
amples of microquasars, accreting black holes or neutron stars in binary systems 
which have associated relativistic jets. A key tool for identifying binary systems 
is their orbital modulation, due to the regular pattern of emission as the orbit is 
traversed. Typically the zero phase of the orbit may be assigned from radio data, 
the time of periastron or some other marker, but is arbitrary. Until 2009, there had 
been no unambiguous detection of y-ray emission from a microquasar system. The 
identification became conclusive with the confirmation of modulation at the orbital 
period. 

We present here the latest results on the HE and VHE y-ray emission from this 
population of sources, with an emphasis on the results and discoveries from the 
Fermi-LAI. 



2 LS I +61 303 

2.1 The original discovery, and further TeV observations: flux, 
spectrum, periodicity 

The very-high energy (VHE) spectrum derived from MAGIC data between 200 
GeV and 4 TeV at orbital phases between 0.4 and 0.7 was fitted by a power law 
function: F y = (2.7 ±0.4 ±0.8) x 10- 12 (£/TeV)- 2 - 6±0 - 2±0 - 2 cm-V'TeV- 1 , with 
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the errors quoted being statistical and systematic, respectively [13]. All further re- 
ports of LS I +61°303 at very-high energies would produce compatible spectra with 
this measurement, within errors. These first MAGIC measurements showed that 
the very-high energy y-ray emission from LS I +61 303 was variable. The maxi- 
mum flux corresponded to about 16% of that of the Crab Nebula, and was detected 
around phase 0.6 with a significance of 8.7c The source is no longer detected near 
periastron and only upper limits were therein imposed. The VERITAS array car- 
ried out independent observations and soon confirmed these results [6]. These early 
MAGIC/VERITAS observations occured at a number of different orbital phases but 
the source was only ever detected at similar phases; this hinted at a periodic nature 
of the emission. 

Further MAGIC observations showed that this emission is indeed periodic [10], 
with the system showing regular outbursts at TeV energies in phases around 0.65 
and no significant signal elsewhere. The periodicity was tested using the Lomb- 
Scargle perodogram both on a set of OFF data and on the source of choice. Whereas 
the former produced no periodicity, the periodogram analysis of the latter did. The 
maximum peak in the periodogram is seen at frequency v = 0.0373d -1 which had 
a post-trial chance probability of 2 x 10~ 7 . Other peaks (although less significant) 
were seen in the signal analysis and not seen in the background sample. When fold- 
ing the data with the frequency obtained, the data did not seem at first sight to be 
described by a simple sinusoidal shape (in part, perhaps due to spacing and quality). 
In fact, fitting a sine function and subtracting it from the power analysis removes 
the peak related to the orbital period but not the others, something that is improved 
when a sine plus a gaussian function is fitted to the data, and subsequently sub- 
tracted from the analysis. Thus, whereas further MAGIC observations succeeded in 
bracketing the periodicity of LS I +61° 303 in TeV y-ray s to be very close to the or- 
bital period (26.8±0.2) days, as well as to other wavelength estimates (e.g., at 1.5a 
from the best estimate, coming from radio of 26.4960±0.0028 days [23]), with the 
error determined by Monte-Carlo simulations, they also reveal a new mystery for 
the flux evolution of LS I +61°303 at these energies. Is there an average (sinusoidal) 
emission along the orbit, plus an extra component associated with the turn on of 
a different process (the Gaussian) that only appears at the TeV maximum phase 
range? 



2.2 TeV and X-ray simultaneous observations 

The first strictly simultaneous observations at TeV and X-rays of LS I +6 1 ° 303 have 
been recently presented by the MAGIC collaboration [28]. Using observations by 
XMM-Newton, Swift, and MAGIC during 60% of an orbit, it was found that there 
is a correlation between the X-ray and TeV emission at the time of the TeV maxi- 
mum. A linear fit to the six MAGIC/XMM-Newton pairs of observations that trace 
the outburst yielded a correlation coefficient of r = 0.97. A linear fit to all ten si- 
multaneous pairs (including Swift data) provides a high correlation coefficient of 
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r = 0.81+021 (which has a probability of about 5 x 10~ 3 to be produced from inde- 
pendent X-ray and TeV fluxes). Due to strong variability of the X-ray emission of 
LS I +61°303, of about 25% in hour scales, other previous campaigns with contem- 
poraneous, but not strictly simultaneous data, as for instance the previous MAGIC 
campaign [11]. 2008, or the VERITAS campaign [5] come short in reaching such a 
result, which, nevertheless, would not suffer from adding further statistics and con- 
firmation. Therefore, the X-ray /TeV correlation found for LS I +61° 303 points to 
emitting processes happening at the same time in both wavelengths, and being the 
result of the same physical population of particles; especially, since at the TeV max- 
imum one would not expect the TeV photons to be subject to significant absorption 
due to gamma-gamma processes. 



2.3 The Fermi results on LS I +61° 303 

With the launch of the Fermi Gamma-ray Space Telescope on 1 1 June 2008 and 
the beginning of operations in its survey mode [14], coverage of the full sky above 
100 MeV is available through its main instrument, the Large Area Telescope (LAT). 
There are two main results to extract out of the Fermi observations of LS I +61° 303 
using the data collected in the first 9 months after launch [1]: 

1 . The unambiguous detection of a source co-located with the position of LS I +6 1 ° 303 
which is seen to be periodic (too) at GeV energies (the first such periodicity ever 
detected) and anti-correlated with the TeV emission. 

2. The discovery of a cutoff in the spectrum of the Fermi (GeV) source. 

LS I +61° 303 is detected at significance of ^70(7 at a position consistent with 
the optical location of the source. The LAT light curve was constructed using aper- 
ture photometry of 100 MeV - 20 GeV events. From a 3.2 hour binned light curve 
the weighted Lomb-Scargle periodogram [20] was calculated. The periodogram ex- 
hibits a clear peak at 26.6 ± 0.5 days (la errors obtained though Monte-Carlo) 
consistent with the known orbital period. The binned LAT light curve folded on 
the nominal orbital period [23] with the usual zero phase at MJD 43,366.2749 [24] 
shows a large modulation amplitude with maximum flux occurring slightly after 
periastron passage (where there is no TeV emission detected, neither by MAGIC 
nor by VERITAS). The folded light curve and the Lomb-Scargle periodogram are 
shown in Figure 1. The overall light curve can be reasonably well fit by a simple 
sine wave. Fitting a sine wave to each of the individual 9 orbits observed finds that 
the best fit amplitude varies between 6.8±0.9 and 2.2±0.9 xl0~ 7 photons cm~ 2 
s _1 [1], which might be indicative of some orbit-to-orbit variability. 

The LAT data was well fit by a power-law plus exponential cutoff returns; the 
chance probability that a simple power-law was incorrectly rejected is 1.1 x 10~ 9 . 
The photon index was found to be T =2.21 ± 0.04 (stat) ± 0.06 (syst); the flux 
above 100 MeV is (0.82 ± 0.03 (stat) ± 0.07 (syst)) x 10~ 6 photons cm~ 2 s" 1 and 
the cutoff energy is 6.3 ± 1.1 (stat) ± 0.4 (syst) GeV. Reduced statistics in a phase- 
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Frequency (days ) Phase 

Fig. 1 The power spectrum (left) and phase folded light curve (right) of LS I +61° 303 at 100 MeV 
- 20 GeV. Left: The weighted Lomb-Scargle periodogram of the LAT light curve. The vertical 
dashed line indicates the known orbital period [23]; the horizontal dashed lines indicate the marked 
significance levels. Right: The phase folded 100 MeV - 20 GeV light curve of LS I +61° 303. The 
dashed lines indicate periastron and apastron of the system. 



bin of 0.1 of the orbit precluded us from finding a statistically significant cutoff in 
each of the bins, and so a simple power law was fitted. No significant variation of 
this photon index with phase was found. 

During part of the Fermi measurements, there was a VERITAS campaign [27] 
observing at higher energies (100 GeV - 30 TeV), and no detection was found even 
when part of the TeV maximum phase region was covered. The lack of strong emis- 
sion from LS I +61° 303 during these observations is at face value somewhat sur- 
prising, especially if one looks at the light curve in the orbit where all are upper 
limits, but does not contradict (yet) previous measurements. The TeV maximum has 
happened around phase 0.6 with excursions down and up in phases of about 0.2 of 
the orbit. The apastron coverage during this campaign was limited, with large data 
gaps, e.g., the coverage between phases 0.72 and 0.82 consists of only 40 minutes. 
Consequently the low statistics of the VERITAS measurements during this cam- 
paign do not allow firm conclusions to be drawn regarding the disappearance of the 
signal or on its deviation from the usual periodic behavior. It is clear, however, that 
any future ground-based LS I +61° 303 campaign will benefit from automatic Fermi 
coverage. At the same time, a strictly periodic behavior would be easy to rule out by 
subsequent MAGIC/VERITAS observations of the TeV maximum, and the problem 
of understanding the origin of the y-ray emission would be even more complex if 
this were so. 
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Fig. 2 High and very-high energy spectral data points of LS I +61°303 (left) and LS 5039 (right). 
Left: Fermi data points are red; MAGIC data are blue (high state phases 0.5-0.7); VERITAS data 
are black (0.5-0.8). Right: The red data points show the spectrum at inferior conjunction (TeV 
maximum, GeV minimum, phases 0.45-0.9; the blue data points show the spectrum at superior 
conjunction (TeV minimum, GeV maximum, phases < 0.45 and > 0.9). Higher energy data come 
from H.E.S.S. Note that the data from the different telescopes are not contemporaneous, though 
they do cover multiple orbital periods. 



3 LS 5039 



3.1 The original discovery, and further TeV observations: flux, 
spectrum, periodicity 

LS 5039 was discovered by the H.E.S.S. experiment to be a y-ray source[8] and 
soon afterwards, periodicity in the TeV y-ray flux, consistent with the known orbital 
period [18], was also reported [7]. Unlike LS I +61°303 H.E.S.S. observations of 
LS 5039 have detected this source at all orbital phases observed. The authors of [7] 
define two broad phase intervals, the inferior conjunction (0.45 < <p < 0.9) and su- 
perior conjunction ( <p > 0.9 and <p < 0.45), for which different (averaged) H.E.S.S. 
spectra were produced. The differential photon energy spectrum (0.2 to 10.0 TeV) 
for inferior conjunction is consistent with a hard power-law with an exponential cut- 
off, where the slope is F = 1 .85 ± 0.06 (stat) ±0.1 (syst) and the exponential cutoff 
energy is at E = 8.7 ± 2.0 TeV (for the fitted function dN/dE ~ E r exp(-£ /E )). 
In contrast, the spectrum for superior conjunction is consistent with a relatively 
steeper (r = 2.53 ±0.07 (stat) ± 0.1 (syst)) pure power-law (0.2 to 10 TeV). The 
statistics collected along the orbit actually allowed to find that the parameters of 
power-law fits to the y-ray data obtained in 0.1 phase binning already displayed 
significant variability. Because of low statistics, and presumably too, of different 
statistics at the higher end of the spectrum in each of these bins, more complicated 
functions such as a power-law with exponential cutoff provided no better fit than 
that of a pure power-law. These results, together with the H.E.S.S. light curve and 
averaged spectra are shown below, where we make a comparison with various mod- 
els. 
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The dataset for the analysis reported by Fermi spanned the period from 4 Au- 
gust 2008, through 29 April 2009, and thus covered multiple orbits of the sys- 
tem [2]. Like LS I +61° 303 this system is also found to be periodic in GeV y-rays 
(3.903±0.005 days) and anti-correlated with the TeV emission. The periodicity is 
thus also consistent with the orbital period. The SED also presents a cutoff at a few 
GeV and is well fitted by a similar function: £~ r exp(— £ /^cutoff), with the photon 
index being T = 1.9± O.l(stat) ± 0.3 (syst); a 100 MeV-300 GeV flux of (4.9 ± 
0.5 (stat) ± 1.8 (syst)) xl0~ 7 photons crrT 2 s~\ and a cutoff energy at 2.1 ± 0.3 
(stat) ±1.1 (syst) GeV. The energy spectra of both LS 5039 and LS I +6T303 are 
shown in Figure 2. Fermi observations were able to further separate into superior 
and inferior parts of the orbit, in phase ranges as defined by H.E.S.S.: GeV observa- 
tions also pinpoint this (GeV / TeV - anti-correlated) flux variability. The fact that 
H.E.S.S. did not find long-term variability makes it reasonable to directly compare 
non-simultaneous H.E.S.S. and Fermi observations. Using the same phase intervals, 
the inferior conjunction slope was found to be T = 2.25 ±0.11 with no exponen- 
tial cutoff energy required to provide a better fit; whereas at superior conjunction, 
a r = 1.91 ±0.16 with a cutoff energy of 1.9 ±0.5 GeV was found. The spectral 
shape in Fermi is softer around periastron (near superior conjunction) and is harder 
around apastron as shown in Figure 3. 



4 Cygnus X-3 

4.1 Historical observations at high energies 

Cyg X-3 is a well known, powerful high-mass X-ray binary with a short orbital 
period of 4.8 hours [15]. The system comprises of a compact object, the nature 
of which is still subject to debate, accreting matter from a Wolf-Rayet compan- 
ion star [36]. The source was discovered in 1966 and has been widely observed 
across the electromagnetic spectrum. It regularly becomes the brightest radio source 
amongst the known binary systems with large flares attributed to its relativistic jets. 
Emission had been detected up to ^300 keV whilst the system shows a complex 
X-ray spectrum which transitions between two main states; 'soft' and 'hard', the 
source is known to flare in radio when entering the 'soft' state with associated rela- 
tivistic plasma ejection events [22]. 

Historically there were reported detections at MeV-PeV energies in the 1970s 
and 1980s. These observations came from 30 MeV - 5 GeV sensitive balloon 
and satellite borne experiments such as SAS-2 and COS-B, from atmospheric 
Cherenkov experiments at 0.1-500 TeV and extensive air shower experiments at 
0.1-10 PeV [15]. However, the detections were typically of low significance and al- 
though some experiments claimed detections, others could not confirm them. These 
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Fig. 3 Fermi average phase folded light curves for LS 5039. Top: 0.1-10 GeV flux variations with 
orbital phase. Bottom: Changes in the hardness ratio, flux(l-100 GeV)/flux(0.1-l GeV), across 
the orbit. 



detections remained doubtful when the next generation of more sensitive ground 
based telescopes failed to confirm the TeV and PeV detection of Cyg X-3 [33, 16] 
and the CGRO-EGRET failed to detect any GeV emission [32]. Hence, claims of HE 
and VHE emission from this microquasar were controversial and highly contested. 
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Fig. 4 Gamma-ray light curves of two of the flaring events observed by Fermi. The LAT fluxes > 
100 MeV (filled circles) are calculated from a likelihood analysis of 4-day long segments of data; 
a power-law with index of 2.7 is assumed in the analysis. The vertical dashed lines indicate the 
periods the LAT detects flaring activity. In the top of each panel the 15 GHz radio flux as measured 
by the AMI and OVRO 40-m radio telescopes is presented; an offset of 0.124 Jy has been removed 
from the OVRO data to compensate for the effect of extended nearby sources that are resolved by 
the AMI interferometer. It is clear that a number of radio flares are detected coincident with the 
y-ray activity detected by the LAT. 
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4.2 The Fermi results on Cygnus X-3 

The analysis of Fermi-LPZY observations of Cyg X-3 [3] spans data taken from 4 
August 2008 to 2 September 2009. The Cygnus region is challenging to analyze 
due to the high levels of diffuse emission combined with the presence of 3 bright 
7-ray pulsars [4]: PSR J2021+4026; PSR J2021+3651; PSR J2032+4127. In fact, 
PSR J2032+4127, lies very close to the location of Cyg X-3 (-30') and contributes 
significantly to the large number of photons detected by the LAT at the location 
of Cyg X-3. The impact of this pulsar was minimized by only using photons which 
arrive during the off -pulse phase. This preserves approximately 80% of the live time 
and produces a significantly cleaner data set. 

The LAT data yield a detection of a point source at the level of 29a at a location 
consistent with that of Cyg X-3. The 4-day binned 7-ray light curve of the source 
flux >100 MeV (sections of which are shown in Figure 4) clearly indicates that the 
source is highly variable. The LAT detects Cyg X-3 during two specific active pe- 
riods: 1 1 October to 20 December 2009 (MJD 54750-54820); 8 June to 2 August 
2009 (MJD 54990-55045). These active phases comprises one or more flares with a 
peak flux of up to —2 x 10~ 6 photons cm _2 s _1 above 100 MeV. The average spec- 
trum over the two active periods is represented by a single power-law model with a 
spectral index r = 2.70 ± 0.05 (stat) ± 0.20 (syst) and an average flux above 100 
MeV of [4.0 ± 0.3 (stat) ± 1.3 (syst)] x 10~ 10 erg cm~ 2 s _1 . Assuming a distance 
of 7 kpc this gives a luminosity, L r ~ 3 x 10 36 erg s _1 , comparable to what was 
reported by SAS-2 [29]. The AGILE mission also reported detecting 7-ray activ- 
ity coincident with the location of Cyg X-3 during these epochs as well as activity 
from 16-17 April 2008 (MJD 54572-54573), prior to the launch of Fermi [35]. The 
AGILE team report a 5.5(7 detection with a comparable average flaring flux to that 
reported by the ferm/-LAT. 

The association of the 7-ray flaring source was definitively confirmed though the 
detection of the 4.8 hour orbital period in a finely binned (1000s) aperture photom- 
etry light curve. The periodicity was found using the Lomb-Scargle periodogram 
and weighting the input flux data points by their relative exposure. No orbital pe- 
riodicity is evident when using the entire data set, however restricting the data set 
to the periods of LAT detected enhanced emission discussed above yields a strong 
detection with a nominal false detection probability of 3.6x10 . Comparing the 
folded 0.1-100 GeV light curve with the 1.5-12 keV X-ray light curve taken from 
RXTE/ASM data shows that both have the same asymmetric shape with a slow rise 
and faster decay. However, the LAT minimum trails the X-ray minimum by 0.3-0.4 
in phase. 

The relationship between the 7-ray and radio emission was confirmed through 
a discrete cross-correlation analysis which indicated a positive correlation between 
the two wavebands with a significance of > 3(7. The lag of the radio light curve to 
the 7-rays is not well constrained and is estimated to be 5 ± 7 days. This specifi- 
cally links the 7-ray activity to periods of relativistic ejection events. The LAT flux 
is compatible with extrapolations of the hard X-ray tail observed above 30 keV, re- 
portedly up to several 100 keV [26], in instances of the 'ultrasoft' state. In this case 
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Fig. 5 Left: The Fermi power spectrum using only the data during the periods of enhanced emis- 
sion (top) and for the entire LAT data set (bottom). The vertical arrows indicate the frequencies of 
the orbital period and the second harmonic. 

Right: The average orbital phase folded light curve of the LAT data during the periods of enhanced 
emission (top) and of the RXTE ASM X-ray data taking during the RXTE mission lifetime. Phase 
zero is set to be at the point of superior conjunction. 



the Fermi emission could be explained by inverse Compton scattering of UV pho- 
tons from the WR star off of high energy electrons. However this scenario requires 
that the emission region is not too close to the accretion disk otherwise the y-ray 
emission will be absorbed through pair production on soft X-ray photons coming 
from the disk. 

If inverse Compton scattering is the dominant y-ray production method then it is 
expected that the peak in the y-rays would correspond to the time of superior con- 
junction when the electrons are seen behind the WR star and hence the energetic 
electrons directed towards the Earth undergo head-on collisions with the stellar UV 
photons. This would correspond to the X-ray minimum assuming the X-ray mod- 
ulation is produced through Compton scattering in the WR wind which is approx- 
imately what is observed in the average folded light curves. A proposed model to 
explain the observed behavior of Cyg X-3 at high energies is that of a jet launched 
around a black hole oriented not too far from the line of sight, which interacts with 
the WR stellar wind to produce a shock at a distance of 1-10 times the orbital separa- 
tion from the system [21]. This shock is the location where electrons are accelerated 
to GeV energies and upscatter star photons. 



5 Cygnus X-l 



Cyg X-l is an exceptionally bright, well studied high-mass X-ray binary and is be- 
lieved to host a stellar mass black hole. It is located at a distance of ~2.2 kpc and 
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comprises of a 09.7 lab donor star in a 5.6 day orbit with the compact object [37]. 
Cyg X-l is another example of a microquasar, it is believed to be principally pow- 
ered by accretion and is seen to transition between 'soft' and 'hard' X-ray spectral 
states. The radio emission is observed to be relatively stable during the 'hard' state 
except for rarely observed radio flares [22] and is undetected in the 'soft' state. 



5.1 Reports of high energy and very high energy emission 

Between June and November 2006 the MAGIC telescope observed Cyg X-l for a 
total of 46.2 hours. The MAGIC team searched the entire data set and found no 
significant persistent VHE y-ray source at the location of Cyg X-l and reported a 
95% flux upper limit in the 150 GeV - 3 TeV band of the order of 1-5% of the 
Crab [12]. Searching for faster time-varying signals within any given night resulted 
in the detection an excess of 4.9a (4.1c after correcting for trials) resulting from 
79 minutes of effective on-time on September 24 2006 UTC 22:17 - 23:41. The 
observed excess is reported to be consistent with a point source at the location of 
Cyg X-l and excludes a nearby radio nebula. The spectrum was fitted by a power 
law function: F y = (2.3 ±0.6) x 10- 12 (£'/TeV)- 3 - 2±0 - 6 cm- 2 s- 1 TeV" 1 , with the 
errors quoted being statistical only; the systematic uncertainty is estimated to be 
35% on the overall flux normalization and 0.2 in the determination of the spectral 
index [12]. 

The TeV flare was coincident with a hard X-ray flare detected by INTEGRAL, 
Swift-BAT and RXTE-ASM. The TeV detection was observed at the rising edge of 
the first hard X-ray peak and the MAGIC non-detection the following night coin- 
cided with the decay of the second hard X-ray peak and hypothesize that the hard 
X-ray and y-rays may be produced in regions linked by the collimated jet [12]; these 
processes would have different physical timescales explaining the shift in time be- 
tween the TeV and X-ray peaks. 

Recently the AGILE satellite has reported the detection of transient y-ray activity 
from Cyg X-l above 100 MeV [34]. They used -315 days of data from July 2007 - 
mid October 2009 giving a net exposure of —13 Ms. Performing a likelihood anal- 
ysis over the entire integrated dataset does not yield a significant detection of Cyg 
X-l and places a 2a upper limit on the y-ray flux (0.1-3 GeV) of 3x 10~ 8 photons 
cm~ 2 s~ 1 . The AGILE team searched for variability on day timescales and reported a 
single flaring episode corresponding to 15 October 2009 UTC 23:13:36- 16 Octo- 
ber 2009 UTC 23:02:24, with a 5.3a pre-trial (4a post-trial) significance level. The 
location of the flare is consistent with the location of Cyg X-l and corresponds to a 
0.1-3 GeV y-ray flux of (2.32 ± 0.66) x 10~ 6 photons cm^s" 1 . At the time of the 
flare Cyg X-l was in the 'hard' X-ray spectral state, however, unlike the instance of 
the MAGIC flaring event there was no simultaneous flare seen in the X-ray (or any 
other) waveband. In March 2010 the AGILE team reported new transient activity 
in the Cygnus region which was compatible with the Cyg X-l location [17]. They 
reported a 5 a detection integrating from 24 March 2010 UTC 02:24 to 25 March 
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2010 UTC 01:01 with a corresponding flux above 100 MeV of (2.0 ± 0.9) x 10~ 6 
photons cm _2 s _1 . 

A search of the Fermi-LXT data for signs of the October 2009 flare yielded no 
significant detection. Extracting data from the exact period reported by AGILE the 
LAT places a 95% upper limit on the y-ray flux (0.1-3 GeV) of ^4x 10~ 7 photons 
cm _2 s _1 (assuming a power-law spectral shape with photon index, T=2.2) [30]. 
Similarly the LAT did not detect any flaring emission during March 2010. 



6 Conclusion 

The first couple of years of Fermi-LAT operations have been highly fruitful in ex- 
panding our understanding of the high energy y-ray emission properties of Galactic 
y-ray binaries. In fact the availability of sensitive telescopes spanning the MeV-TeV 
energy range in the past decade have opened up this entire field of study and the LAT 
is the latest instrument to add to and expand this field of study. 

Gamma-ray binaries have now been definitively detected in the GeV energy 
range; both LS I +61° 303 and LS 5039 are strong detections which have been iden- 
tified through their spatial location and the clear signs of orbital modulation in the 
LAT observed flux. The orbital modulation of the GeV flux was found to be anti- 
correlated with the modulation observed at TeV energies with the GeV flux peaking 
around periastron in both cases. The identification of an exponential cut-off power- 
law shape to the spectrum of both sources was not expected. It has been noted that 
the spectral shape observed in both y-ray binaries is reminiscent of that seen in the 
large family of y-ray pulsars now detected by the LAT [2] (i.e. a hard power-law 
spectrum with an exponential cutoff at ^2.5 GeV). It has been proposed that this 
may indicate that the emission in the Fermi energy range from both LS I +61° 303 
and LS 5039 may therefore be magnetospheric in origin, however, this raises a new 
problem as it is not clear how such emission would be modulated by the orbital pe- 
riod of the binary. A clear indication of magnetospheric emission would come from 
the detection of pulsations in the Fermi emission of these sources whilst simultane- 
ously confirming for the first time the nature of the compact object. 

PSR B 1259—63 is the only known millisecond pulsar in a binary system with a 
main-sequence star. It has been detected at TeV energies by H.E.S.S. [9] during peri- 
astron passages of its 3.4 year orbital period. Fermi has yet to report any detection of 
pulsed magnetospheric emission from this source however, the upcoming periastron 
passage in December 2010 will offer the first opportunity for Fermi to observe the 
expected GeV emission resulting from the pulsar wind interacting with the stellar 
wind of the companion. Observations of this periastron passage will provide a key 
dataset to compare and contrast with those of the persistent y-ray binaries in which 
the nature of the compact object is still unknown. 

The detection of transient y-ray activity from Cyg X-3 by both Fermi and AG- 
ILE during multiple epochs of soft-state transitions has confirmed this microquasar 
as a source of high energy emission. Furthermore, it has been seen that the y-ray 
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emission is correlated with radio flare events. The clear signature of orbital flux 
modulation detected by the LAT removes any doubt that Cyg X-3 is the source of 
the emission whilst also suggesting that the emitting region is likely bounded by the 
orbital separation of the system. 

The recent reports by AGILE of y-ray activity from another microquasar, Cyg 
X-l, have not been confirmed by Fermi. The potential to discover flaring events 
from this source and other microquasars remains a real possibility through future 
observations with the LAT. Indeed with the upcoming periastron passage of PSR 
B 1259— 63, the detection of a symbiotic binary undergoing a nova outburst [19] and 
the ever present possibility of detecting a new y-ray binary the future is very bright 
for Fermi-LKT observations of binary systems. 
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